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Abstract 



We explore possible CP violating effects, coming from the Dirac phase of 
(— | \ the Maki-Nakagawa-Sakata-Pontecorvo matrix, on the neutrino degeneracy 

Qh' parameters, at the epoch of Big-Bang nucleosynthesis. We first demonstrate 

O ■ the conditions under which such effects can arise. In particular it requires 

that the initial muon and tau neutrino degeneracy parameters differ. Then 
we solve numerically the kinetic equations for the three flavour neutrino 
density matrix with the goal of quantifying the impact of the Dirac phase 
on £ Ve . The calculations include the vacuum term, the coupling to matter, 
the vv interaction and the collisions. Effects on £ Ue up to almost 1.% and 
on Y p of about 0.1% are found, depending on the initial conditions. 

Keywords: Physics of the early Universe, Neutrino physics, CP violation 
g ; PACS: 14.60.Pq, ll.30.Er, 26.35,+c 
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1. Introduction 

^ ■ One of the major open questions in modern cosmology is the origin of the 

matter-antimatter asymmetry in our Universe. The bayon asymmetry is 
nowadays known to be rjB = (ns— ng)/n 7 = 6.14x 10 _10 (1.00±0.04) thanks 
to the measurement of the CMB anisotropies by WMAP [lj . Sphalerons ef- 
fects in baryogenesis and leptogenesis scenarios [3] can equilibrate cosmic 
lepton and baryon asymmetries at the same level. Since the lepton asym- 
metry is only possible in the neutrino sector because of charge conservation, 
the observation of a non-zero neutrino degeneracy parameter £ can furnish 
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important information to our understanding of the matter-antimatter asym- 
metry in the Universe. 

In analogy with rjs related to the baryon asymmetry, the total neutrino 
asymmetry L u = L Ue + L„ + L Vr can be quantified by the neutrino chemi- 
cal potentials fi Ua (a = e, /i, r) or, equivalently, the degeneracy parameters 

j _ n Vg ~ n Vg _ 71-2 ( r ^Vg \ (t , ^Vg \ I -i\ 

Lva ~ n 7 " 12C(3) \ Tj V^ q vr2 J W 

where n Va {n Va ) are the neutrinos (anti-neutrinos) occupation numbers and 
C(3) ~ 1.202. Non-zero electron, muon and tau neutrino degeneracy pa- 
rameters influence the abundance of light elements produced in Big-Bang 
Nucleosynthesis (BBN) in two aspects. While all flavours influence the ex- 
pansion rate of the Universe, by modifying the effective number of degrees 
of freedom, only £ ue impacts the neutron/proton ratio, a key parameter for 
the 4 He abundance. Indeed 4 He, among all the light elements formed during 
BBN, is the most sensitive one to the neutrino degeneracy parameters. Ex- 
tensive work has been performed to extract information on the relic lepton 
asymmetries either from Big-Bang Nucleosynthesis, as in e.g. @)0,@)S,Q)Sl 
or from the cosmic microwave background and large scale anisotropics, like 
in @. 

Major advances have been performed in neutrino physics in the last ten 
years. The change in neutrino flavour content due to oscillations is at present 
a well established phenomenon. This implies that the neutrino flavour basis 
is related to the mass basis 



where the unitary Maki-Nakagawa-Sakata-Pontecorvo (MNSP) matrix can 
be written as a product of three matrices U = T23T13T12 

U= ^0 C 23 S23 I I 1 || -8 12 C 12 | , (3) 

with Cij = cosOij (sij = sinOij) and 9x2, 623 and #13 the three neutrino mixing 
angles. These oscillation parameters have been well determined, except 
for the third neutrino mixing angle #13 and a possible Dirac CP violating 
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phasqj. If #i3 is close to the Chooz limit, i.e. sin 2 2#i3 < 0.02, reactor 
experiments (Double-Chooz, RENO and Daya-Bay) should soon measure 
this angle floT ]. The two squared mass difference^! have been measured with 
good precision Since the sign of Am^ has not been determined yet, two 
mass hierarchies are possible: inverted (Am^ > 0) or normal (Am^ < 0). 
This is known as the mass hierarchy problem. The absolute neutrino mass 
scale is also still unknown, since neutrino oscillations are only sensitive to 
mass squared differences. The KATRIN experiment will soon reach the sub- 
eV sensitivity 12] while important indirect limits on the sum of the neutrino 
masses are obtained using CMB and LSS data (see e.g. 012 EH). Indeed, 
so far, only indirect effects of cosmological neutrinos have been observed. 
Their detection represents one of the major future challenges. An interesting 
possibility has been proposed recently in [iff], namely to exploit the capture 
on radioactive nuclei. This idea has been further investigated in 171. Il8j|. 

The observation of CP violation in the neutrino sector is a key open 
question. The breaking of the CP symmetry can arise from the presence of 
a non-zero Dirac 5 phase that renders the U matrix complex Eq.©. Long- 
term expensive accelerator complex (super-beams, beta-beams or neutrino 
factories) might be required to tackle this issue [191 ] ■ Therefore it is impor- 
tant to explore complementary avenues and search for indirect effects. For 
example, recently we have explored possible CP violation effects in core- 
collapse supernovae [2(|. We have shown that they can arise e.g. if i/„ 
and v T experience a different refractive index in the medium - due to loop 
corrections and/or physics beyond the standard model. These results have 
been extended in presence of the neutrino-neutrino interaction in 2l| . Note 
that important developments are currently ongoing in the study of neu- 
trino propagation in dense media due to temporally evolving density profiles 



24l. |25|. l26l. |27I]. and the neutrino-neutrino interaction, which intro 



duces collective phenomena (see e.g. [28|, |29j, |30j, [3l|, |32|]). The importance 



of the latter contribution has been first pointed out in the early Universe 



context 



34, 



Several calculations have been performed of the neutrino degeneracy evo- 
lution at the BBN epoch including neutrino oscillations 
In 3j] it is shown that neutrino oscillations tend to equilibrate the elec- 



1 Note that Majorana phases can also be present. They can influence the neutrinoless 
double-beta decay half-lives, while neutrino oscillations are not affected by such phases. 
For this reason they will not be considered here. 

2 Although the existence of sterile neutrinos is an attractive possibility, here we consider 
three active neutrino families, in agreement with the ensemble of experimental data. 
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tron, muon and tau neutrino degeneracies. However how much flavour equi- 
libration really holds is still unclear [33]. Only in case of flavour equili- 
bration the constraints on £ Ve coming from the abundance of primordial 
Helium-4 can be translated to the other flavours. In such a case the limit 
of —0.044 < £ < 0.070 for all flavours 0], if the conservative Olive and 
Skillman analysis is used, with an uncertainty of the order of 5% on Y p [38| . 



Other detailed analysis of the 4 He fraction exist 39|, 140J. If the systematic 
uncertainties, inherent to the helium abundance measurements, are better 
understood a precision as low as 10 -3 might be reached. Besides, future 
studies of gravitational lensing distortions on both the temperature and the 
CMB polarization might reach sensitivities, close to the BBN ones, on the 



helium fraction, and even at the level of 5. 10~ 3 41]. 

In this paper we explore possible CP violating effects, coming from the 
Dirac phase, on the neutrino degeneracy parameters, at the time of Big- 
Bang Nucleosynthesis. First we demonstrate analytically the conditions 
under which there can be such effects. Then we determine numerically 
the neutrino degeneracies evolution including, for the first time, a non-zero 
Dirac phase. We solve the equations for the three flavour density matrix 
taking into account the vacuum oscillations, the coupling to the plasma, the 
neutrino-neutrino interaction and the collisions, using a damping approxi- 
mation. 

The manuscript is structured as follows. The theoretical formalism is 
shown in section 2. Section 3 presents our analytical results that define 
the conditions to have CP effects on the neutrino degeneracy parameters. 
Section 4 illustrates the numerical results and the potential modifications 
introduced by the Dirac phase as well as a discussion of the possible impli- 
cations on the 4 He fraction. Section 5 is the conclusion. 

2. The neutrino evolution equations 

The neutrino evolution including oscillations can be determined by using 
the density matrix: 

Pu(P,t) = I Pv, e Puuu Pvut I ( 4 ) 







Pver 






PV^T 


Pv Te . 


Pv Ti i 


Pv TT 



in three flavours. Each neutrino state is characterized by the momentum 
p and the time t. The transposed of p u (p,t), p u (p,t), is taken to describe 
anti-neutrinos. In an expanding universe, the equations of motion are [291 ]: 

i(d t - Hpd p )p p = [H tot , Pp] + C(pp), (5) 
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where the explicit dependence on t is not shown for simplicity, the subscript 
p refers to the momentum dependence, and Htot is the total Hamiltonian 
describing neutrino propagation in the medium. As long as the expansion 
rate of the Universe is smaller than the collision rate among the relativistic 
species, collisions play an important role and drive the system towards equi- 
librium. Such contributions, proportional to G\ are included here through 
the collision term C(p). The cosmic expansion is taken into account through 
the Hpdp contribution, with H = a(t)/a(t). a(t) is the scale factor that is 
normalized such as a ~ 1/T at high temperatures or early times. The Hubble 
constant H is determined through the Friedmann equation H = ^J&irGp/?> 
with G being the gravitational constant and p the total energy density of 
the relativistic particles. By using co-moving variables x = ma, y = pa (m 
is an arbitrary mass scale that we take equal to 1 MeV) Eq.([5]) becomes 
adimensional: 

lHxd x Py = [H tot ,Py] + C(Py) (6) 

The total Hamiltonian describing neutrino propagation involves three 
contributions 



UM 2 W 8V2G F p 



Htot = H vac + H mat + H vv = — — 2 — E + V2Gf(p — p), (7) 



where Gf is the Fermi constant and mw the W boson mass. Anti-neutrino 
evolution is described by the same Eq.([7J) for p u (p,t), but with a minus 
sign for H vac . The first term is the vacuum oscillation contribution with 
M 2 = diag{m\,m\,m\) and U the MNSP matrix Eq.©. The second contri- 
bution is proportional to the energy densities E of charged leptons (electrons, 
positrons and muons) in the plasma and corresponds to the refractive effects 
of the medium that neutrinos experience [29(. Note that the background 
potential arising due to asymmetries in charged leptons is negligible in com- 
parison with the other terms in the Hamiltonian [j^ ]. The \f2Gp(p — p) 
contribution represents the neutrino-neutrino interactions and is meant to 
be integrated over the neutrino momenta. This non-linear term is respon- 



sible for synchronizing the neutrino ensemble, as discussed in [35|, |34|, l33f j . 
similarly to what occurs in core-collapse supernovae 28|, 3(], 21]. Concerning 
the neutrino scattering with e^, p^ or among themselves, in principle one 
should consider the exact collision integral I Va 4^, 37] including all relevant 
two-body weak reactions of the type u a (l) + 2 — > 3 + 4. Here we follow 
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34, 36, 37 1 and use a damping prescription of the forrd_| 

C(Py,a^) = —D a /3Py,aP (8) 
C{py, aa ) = D aa (f(y,^ a ) — Py^, 



\aa J 



with a, j3 = e, p, r and £ a being the equilibrium solution. The coefficients 
are fixed at D a p = 2(Asin 4 6\y — 2sin 2 9w + 2)Fo (9w being the Weinberg 
angle) for a = e and j3 = e, p or r, while for all other cases we take D a p or 
D aa = 2(2sin A 6 w + 6)F with F as in 0]. 

We consider the plasma to be in thermaQ but not chemical equilib- 
rium. Therefore, before making the density matrix evolve, we consider the 
neutrino occupation numbers given by Fermi-Dirac distributions f(y, ^J, 
characterized by the temperature T and the chemical potentials 4j : 

/ f(y,U) o o \ 

p u (y,t = 0)=\ f(y,^) (9) 
\ f(y,U) J 

Before showing the possible impact of the CP phase 5 on ^, we now try 
to get an analytical insight on the possible sources for these effects. 

3. Conditions for CP effects on £ : Analytical results 

Let us now demonstrate under which conditions there can be CP violation 
effects coming from the Dirac phase at the BBN epoch. To this end, we fol- 



low partly the procedure established in Refs.[20] and 2l| within the context 
of core-collapse supernovae. 

For our purpose it is convenient to work in the T23 basis, as shown in 
Ref |2o| . but applied to the density matrix Eq.([3]) 

Pv,y = T 23Pv,y T w (10) 

We also define the useful quantity p y ^s = S^p y S where the 5 dependence is 
contained in the unitary diagonal matrix 



5= 1 (111 




3 Note that in [3y, [33 the damping prescription is used only for the off-diagonal con- 
tributions. 

4 Corrections to the Fermi-Dirac distributions have been calculated to be very small 
5 Opposite chemical potentials are considered for v and v. 
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Since the MNSP matrix Eq.Q can be written as U = T2^ST^T\2, we 
multiply Eq.© by T^ 3 (T23) on the left (right) and get 



iHxd x p Vj s 

with 



Htot,Py,S +C(p y , S ) (12) 



- m = X^M^L _ ZgGpM + _ ^ (13) 

where E 1 = Trjg E T23 and E 1 = diag{E ee ,E IJiljL ,Qi). The quantities _E ee and 
are the energy densities associated with the electrons, positrons and 
p + , fi~ respectively. 

We now show the conditions under which : i) the initial conditions for p$ 
and p are the same, ii) the evolution equations for ps Eq. (|12p are the same 
as for p Eq.([6]). If both i) and ii) are fullfllled then using time discretization 
by mathematical induction one can show p$ is equal to p at all times, and 
therefore the density matrix does not depend on 5 at any time. On the 
contrary, there can be CP violating effects on 5 and on £. 

Let us consider the initial conditions for ps and for each of the quantities 
on the r.h.s. of Eq. (|12p to identify the conditions under which ps = p(t = 0), 
namely when condition i) is satisfied. At t = the density matrix is 

/ f(v,U) 

Py,S=\ c| 3 /(y,^) + s| 3 /(y,^ T ) c 23 s 2 3(/(y, &-) - f(y, 6/„))e* 

\ c 23 s 2 3(/(y,U) - f(y,^))e~ iS 4sf(v,U) + 4»/(i/>&„) 

(14) 

This implies that, if the initial muon and tau neutrino degeneracy parame- 
ters differ, this engenders a dependence on 5 of the density matrix. Concern- 
ing the different terms on the r.h.s. of Eq. (|12|) . it is first shown in Ref.(2o| 
that the vacuum contribution to the Hamiltonian satisfies the factorization 
H va c{3 = 0) = H vac (6)S and therefore has no 5 dependence at any time. 
The matter related term in Eq. (|12j) is given by: 



/ E ee \ 

SES\t = 0)=l -4n E w . c 2 3S 23 E^e tS (15) 
\ c 23 S23E fMfl e- i5 -<%3 E w J 

One can see that, if the presence of muons and anti-muons in the relativistic 
plasma is not neglected at this epoch of the Universe evolution, then this will 
introduce a source of CP-violation. Note that it has been explicitly shown 
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in Ref. 44] that a difference in the and v T refractive indeces induces a 
dependence on 5 in the v e channel. However, the more the temperature 
goes down, the less the E„„ term will be important. Concerning the vv 
interaction contribution at initial time, one has before integrating over the 
neutrino momenta : 

/ L VK \ 

Ps-~Ps = n i\ c 23 L i> + S 23 L ^ c 23 s 2 3(^ - L Ur )e lS J (16) 
\ c 23 s 2 3(^ M -L Ur )e~ iS sl 3 L Ufi + (% 3 L Vt J 

where L Ui is the i flavour lepton asymmetry and n 7 the photon number 
density. Any lepton flavour asymmetry between muon and tau neutrinos 
introduces a CP dependence of the neutrino-neutrino interaction Hamilto- 
nian. Finally, one has that the collision term C(p(p, t)) has no 5 dependence 
since p vm = f(y,£,i) and the p v . Vj = for i ^ j making such a term equal 
to zero in any basis initially. 

Let us now show that condition ii) holds. As shown in [2Qj the mat- 
ter term is equal in both equations if and only if one can neglect the E^ 
contribution. Concerning the neutrino-neutrino interaction term, since the 
corresponding Hamiltonian has a linear dependence in the density matrix 
S^H UV {5)S = H vu {5 = 0) at all times, as demonstrated by mathematical 



induction in Ref. 21], using the Liouville-Von Neumann equation, in the 
supernova context. Here this requires the muon and tau neutrino asymme- 
tries being equal (L Ufi = L Ut ) at initial time. Finally the collision term, 
here treated in the damping approximation, is also a linear function of the 
density matrix, therefore its dependence is the same if ps and p are used. 

In conclusion, if is negligeable and L v = L Ur at t = 0, then the 
Hamiltonian that governs the evolution of the density matrix with and with- 
out a dependence of the Dirac phase is the same. As a consequence, the 
evolution of ps and p is the same. This implies that, under such conditions, 
PvevA$) = Pu e uA& = 0), and therefore £ Ve {5) = £v e {?> = 0) at all times. 

Summarizing, we have demonstrated that possible CP violating effects 
can arise from the Dirac phase 5 if, initially, there is a difference between the 
muon and tau neutrino occupation numbers and/or degeneracy parameters. 
If the (anti)muon contribution to the energy density is non-negligible this 
can also, in principle, engender CP violating effects. 

4. Impact of the CP phase on ^i/ e ' Numerical results 

The goal here is to quantify possible CP effects on the neutrino degenera- 
cies just before Big-Bang nucleosynthesis, having in mind that this can e.g. 



S 



potentially impact the Helium-4 fraction. The numerical results we present 
are obtained by solving Eqs.([6]l8|) for the density matrix and the initial con- 
ditions Eq.©. The neutrino degeneracies are then obtained using Eq.([T|). 
We will give the variations on £„ e induced by S since this is the relevant 
quantity to quantify the effect on the helium-4 fraction. 

In our calculations, the oscillation parameters are fixed at the values 
Am 2 12 = 8x 10" 5 eV 2 , sm 2 28 12 = 0.83 and Am| 3 = 3x l(T 3 eV 2 , sin 2 20 2 3 = 1 
for the solar and atmospheric differences of the squared mass differences and 
mixings, respectively 111 ] . For the third still unknown neutrino mixing angle 
^13, we have taken either a large value close to the Chooz limit, namely 
sin 2 20i 3 = 0.19 at 90 % C.L., or a small value of sin 2 26>i 3 = 3 x 10 -4 . 

Figures HH] present the evolution of the neutrino degeneracy parameters 
as a function of the temperature. Figure [T] shows the results obtained for 
£ with 5 = 0° and 5 = 180°, without the inclusion of the neutrino- neutrino 
contribution to the total Hamiltonian Eq. <113 j) . Only the off-diagonal contri- 
butions to the collision term are included. One can see that the effect of S 
are very small, the curves for a non-zero 5 value being indistiguishable from 
those for 5 = 0°, as we have been verifying for different initial conditions. 
For example, if at initial time £„ e = —0.3, £„ = 0.3 and £ Ut = 0, the varia- 
tion induced by the phase is A£j, e = 6, e (<5) — ^ e (S = 0) = 4 x 10~ 6 , while 
for £„ e = -0.5, t, Ufl = 0.5 and £„ T = 0, we obtain A£„ e = 10~ 5 . 

Figures [2] and [3]|4] show the CP effect on £ without and with the neutrino- 
neutrino contribution^], respectively. Here both off-diagonal and diagonal 
collision terms are includeqj. Results are shown with 5 = 0° and 6 = 180°. 
In Figure [2] the effect of the CP phase on ^ Ue is at the level of 10~ 3 . The 
degeneracy parameters are compatible with the bounds valid if total flavour 
equilibration is assumed. In Figj3]the effect is A£j, e = 6 x 10~ 3 . Note that, 
for such initial conditions, while the value of £„ e at decoupling is compatible 
with the present bound from BBN; the other degeneracy parameters are 
compatible with the bounds valid without assuming total flavour equilibra- 
tion. Figure H] presents A£„ e for different initial conditions, showing that 
the variations induced by 5 increase, if £ Ve at the temperature of neutrino 
decoupling is larger. 

In all the calculations performed we have found that the modifications 
produced by the CP effects increase for a maximal phase and if the third 



6 The calculations including the neutrino-neutrino contribution have been performed 
following the procedure used in [§3, EH ■ 

7 We use the approximation that the equilibrium £ are kept at the initial value [34| . 
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Figure 1: Neutrino degeneracy parameters as a function of the temperature for 5 = 0° 
and 5 = 180° , at the BBN epoch. The initial conditions are set at £„ e = 0.1, ^ Vt = —0.1 
and £„ = 0. Results obtained solving Eqs. (|6I8[I and the initial conditions Eq.Q but 
without including the neutrino-neutrino interaction and the diagonal collision terms. The 
third neutrino mixing angle is taken at the Chooz limit. Here the curves corresponding 
to a non-zero Dirac phase are indistinguishable from those with a zero value. 



neutrino mixing angle is large. Note that the CP effects arising from the 
presence of a non-zero have been found to be completely negligeable. 
The inclusion of the non-linear vv contribution, in general, reduces the phase 
impact. This can be qualitatively understood. In fact such term synchro- 
nizes the neutrino ensemble and freezes neutrino flavour conversion. A simi- 
lar reduction of CP effects on the neutrino oscillation probabilities has been 
found in the context of core-collapse supernovae, when the neutrino cou- 
pling to neutrinos is included (see Figure 3 of Ref . [2JJ] ) . Finally, our results 
show that the effect depends on the diagonal collision terms. Clearly, a def- 
inite conclusion on the quantitative effects needs the inclusion of the exact 
collision integrals. This will be the object of future investigations. 

As far as the impact on the on 4 He fraction is concerned, an estimate of 
the CP effects can be obtained using the simple relation AY p = —0.2A^ Ve 
[itH - In fact, a modification of the order of A£ e of a several 1CP 3 , as we 
have found for some initial conditions, modifies Y p at most by about 10 -3 . 
This is within the uncertainty from BBN observations, although the present 
uncertainty on the 4 He fraction might be narrowed down at the level of 
Y p < 0.005 in the future [II]. 



10 



30 20 10 5 2 

Temperature (MeV) 



Figure 2: CP effects on the £ as a function of the temperature. The initial conditions 
here are £„ e = £„ T = 0. and £ v = —0.1. The calculations include collision terms but do 
not include the neutrino-neutrino interaction. The mixing angle #13 is taken at the Chooz 
limit. 

5. Conclusions 

We have explored the impact of the Dirac CP phase of the MNSP matrix on 
the neutrino degeneracy parameters at the Big-Bang nucleosynthesis epoch. 
First we have established analytically the conditions under which there can 
be possible CP effects coming from this phase and shown, in particular, that 
these are present if there is a difference between the initial muon and tau 
neutrino degeneracy parameters. To quantify such effects we have numeri- 
cally solved the evolution equation for the density matrix in three flavours, 
including mixings in vacuum, coupling to matter, the vv interaction and 
collisions (in the damping approximation). We have found that, depending 
on the initial conditions for the neutrino degeneracy parameters, which are 
an unknown, modifications up to almost 1.% and 0.1% might be present on 
£„ e and Y p respectively, when the vv interaction and the collision terms are 
included. 
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Figure 3: Neutrino degeneracy parameters, as a function of the temperature, when the 
initial conditions are taken equal to £„ e = = and £ Vt = 0.5. The results correspond 
to £„ c for S = 180° and <5 = 0° (lower lines), to £„ T for 5 = 180° (dashed) and 5 = 0° (dot- 
dot-dashed) and to for 5 = 180° (dotted) and 5 = 0° (dot-dashed). The calculations 
include the vacuum oscillation and matter term, the neutrino-neutrino interaction and the 
collisions. 
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